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Unlike general relativity, scalar-tensor theories of gravity predict scalar gravitational waves even 
from a spherically symmetric gravitational collapse. We solve numerically the generation and propa- 
gation of the scalar gravitational wave from a spherically symmetric and homogeneous dust collapse 
under the approximation that we can neglect the back reaction of the scalar wave on the space-time, 
and examine how the amplitude, characteristic frequency, and wave form of the observed scalar 
gravitational wave depend on the initial radius and mass of the dust and parameters contained in 
the theory. In the Brans-Dicke theory, through the observation of the scalar gravitational wave, it 
is possible to determine the initial radius and mass and a parameter contained in the theory. In the 
] scalar-tensor theories, it would be possible to get the information the first derivative of the coupling 

function contained in the theory because the wave form of the scalar gravitational wave greatly 
depends on it. 
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I. INTRODUCTION 



Until now general relativity has passed all the test experiments in the solar system and pulsar-timing tests with high 
accuracy. However, general relativity is not the only gravitational theory that passes such tests. There are alternative 
theories of gravity that pass the present weak-field gravitational tests. In the strong-field regime, these theories may 
deviate very widely from general relativity. Therefore experiments in strong-field regimes are necessary in order to 
determine the correct theory of gravity. 

In the 1960s scalar-tensor theories of gravity were actively studied as alternative theories of gravity. Recently they 
have been revived. The reasons are that these theories may play an important role in the hyperextended inflation 
model that these theories arise naturally as the low-energy limit of string theory |^,|| or other unified theories, 
and that the laser interferometric gravitational wave observatories (LIGO project will be in practical use in the 
next decade and other projects (VIRGO §], GEO 0, and TAMA g) of gravitational wave observations are being 
carried forward so that high-accuracy tests of the scalar-tensor theories may be expected ||[l0) . 

In the scalar-tensor theories of gravity we consider scalar fields as well as a metric tensor as fields that are related 
to gravity. The scalar-tensor theories must agree with general relativity within an accuracy of at least ~ 0.1% in 
the weak- field region from the results of solar-system experiments and pulsar tests . But a class of scalar-tensor 
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theories that satisfies the above constraints may predict order-of-magnitude deviations from general relativity in the 
regime of strong gravitational field One example of such phenomena in the strong-field regime may be the 

gravitational collapse of a star. The Brans-Dicke theory (ll| is the well-known example of scalar-tensor theories. 
Collapse in the Brans-Dicke theory has been researched in comparison with that in general relativity. Matsuda and 
Nariai studied a spherically symmetric collapse of ideal gas by numerical calculations of relativistic hydrodynamics 
in the Brans-Dicke theory. Shibata, Nakao, and Nakamura |1C}] solved numerically a spherically symmetric collapse of 
inhomogeneous dust and examined the character of the gravitational collapse and the wave form of the radiated scalar 
gravitational wave. They found that the final state of the collapse is the Schwarzschild black hole, which is consistent 
with the theorem of Hawking |l5| that a stationary black hole in the Brans-Dicke theory is identical with that of 
general relativity. Scheel, Shapiro, and Teukolsky [[l6) found similar results as Shibata, Nakao, and Nakamura |fuj| . 

Since in general relativity no gravitational wave is radiated from a spherically symmetric collapse, we cannot obtain 
any information of collapse from observations of gravitational waves if realistic gravitational collapses to black hole 
are almost spherically symmetric. However in scalar-tensor theories scalar gravitational waves are radiated even 
from spherically symmetric collapse and these scalar gravitational waves can be detected by laser interferometric 
gravitational wave detectors ju}. It is expected that these scalar gravitational waves reflect directly not only the 
stellar initial radius, mass and density but also the theory of gravity. 

In this paper we take a spherically symmetric, homogeneous dust collapse as a model of the collapse of the stellar 
core and examine the scalar gravitational wave radiation and its wave form. As a first step to understand the nature 
of the scalar gravitational waves in general, we neglect the back reaction of the scalar field on the space-time in this 
paper. The method for numerical calculations is similar to that of Cunningham, Price, and Moncrief [Ti| . We use 
null coordinates in order to see the propagation of the scalar gravitational wave. 

This paper is organized as follows. In Sec. || we review scalar-tensor theories of gravity and scalar gravitational 
waves. In Sec. [II we explicitly present our approximation and show the basic equations used in our numerical 
calculation . In Sec. [V we show numerical results and their implications. Sec. [V] will be devoted to discussions. 

We use units of c = 1. The Greek indices fi, v, ■ ■ ■ run over 0,1,2,3. A comma represents a partial derivative. We 
follow Misner, Thorne, and Wheeler's [fl8f sign conventions of the metric tensor, Riemann tensor, and Einstein tensor. 



II. SCALAR-TENSOR THEORY 



A. Field Equations 

We consider a class of scalar-tensor theories of gravity in which gravitation is mediated by a mass-coupled, long-range 
scalar field in addition to space-time curvature. The action is given by |]llf 

/= I(b/ ^^-^^M^x + I^^g^}, (2.1) 

where g^ v is the metric tensor, R is the scalar curvature, and <f> is the scalar field that couples to gravity. <j) approaches 
the cosmological value 0o m the asymptotic region. is a dimensionless arbitrary function of </>. ^ m represents 

matter fields, and I m is the action of them. 
Varying the action by g^ and </> we have 



+^(V / ,V„^-ff / ,„Cty), (2.2) 

u * = i^ ) ^ T -%^^)> (2 - 3) 

where V v and □ are a covariant derivative and the d'Alembertian of g^ respectively. The equations of motion for 
matter fields are given by 

V„T^ = 0, (2.4) 

where the energy-momentum tensor T^ v is defined by 

rpjxv _ 2 5I m [^m,g^v\ ^2 
V~5 S 9n» 
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The weak equivalence principle is valid in this theory: A test particle which has negligible self-gravitational energy 
moves on a geodesic in this frame, which is often called the Brans-Dicke frame or the Jordan-Fierz frame. However, a 
self-gravitating body moves on a different trajectory because of the coupling of the scalar field to the body's self-gravity 
and therefore the gravitational weak equivalence principle is not valid (Nordtvedt effect) Jll|,|l^] . 

In the asymptotic region, the gravitational "constant" Go, which is measured in time dilation experiment or by 
observation of the Keplerian motion of a test particle, is given by PQ| 



L _ 1 4 + 2io((t)o) 
3 + 2u){(j) Q )' 



Taking the units of Go = 1, we have 



MAO m 



2w(0o) 



If u>((j>) is a constant, the theory reduces to the Brans-Dicke theory. In the limit of uj — > oo, the theory approaches 
general relativity, although it does not mean that every solution of the theory approaches some solution of general 
relativity in this limit [n6|H] ■ 

We can take another representation of the theory by the conformal transformation from the "physical" Brans-Dicke 
frame to the "unphysical" frame in which the Einstein-Hilbert action is separated from the scalar field action p2[ . 
This new frame is often called the "Einstein" frame. In this frame, the action is [E3fl 
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167TG* 



^gl{R» ~ 2g^ u (p^(p >u )d i x + I m [^ m ,A 2 (ip)g^], (2. 



where G* is the "bare" gravitational constant and p is the scalar field. The metric tensor <?» M „ in the Einstein frame 
is related to the metric tensor g^ v in the Brans-Dicke frame by the conformal transformation as 

g llv = A 2 (ip)g* llv . (2.9) 
i?„ is the scalar curvature of <?* M „. tp and A(ip) are determined by 

G*A 2 (p) = I (2.10) 



where a(p) is defined by 



The field equations are given by 



3 + 2cj(0)' 



_ d\nA(p) 

= • (2-12) 

dip 



G*^v = R*nv - -^R*g*^v = 87rG*T„ M ^ + 2 (u>^, v - -^g*ij,ug* 'p,a ( p,0j , (2.13) 

D m tp = -4nG*a(ip)T*, (2-14) 
with the equations of motion given by 

V.„T^ = a(¥>)T*V*^, (2.15) 

where V w , □» and i?» p „ are a covariant derivative , the d'Alembertian and the Ricci tensor of g*^, respectively. 
T* u is defined by 

T *» v _ 2 SI m [^ m ,A 2 {p)g^ l/ ] ^ 

V~g* s g*n» 

and is related to the energy-momentum tensor in the physical Brans-Dicke frame by 
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T^ v = A 2 (ip)T^. (2.17) 

ip approaches the cosmological value <po in the asymptotic region. Hereafter we take the units G* = 1 . 

Generally, a(ip) is an arbitrary function that characterizes the theory. If a(ip) is constant, the theory reduces to 
the Brans-Dicke theory. In the limit a — * , the theory approaches general relativity. We denote 

ao = «W, (2.18) 
da((p) 



1% 



dip 



(2.19) 



(p=ip 

In the post-Newtonian approximation, the PPN parameters ^Edd and 0Edd ar e expressed as 

2al 



l-lEdd=r r r L 5> ( 2 - 2 °) 



l + «0 
2(1 + a%y 



P Ed d-l=^%- 2 . (2.21) 



Results of solar-svstem test experiments (the time-delay and deflection of light) and pulsar-timing tests constrain 
|1 - lEdd\ as Jll|Q 

|l-7B«w| <2x 1(T 3 . (2.22) 

This inequality is rewritten in terms of u> and ao as 

lj > 500, a 2 < 10~ 3 . (2.23) 

The deviation of ^Edd from unity implies the breakdown of the gravitational weak equivalence principle. Bounds on 
this effect from the results of the lunar-laser-ranging experiment constrain \0Edd — 1| as [ p5[ 

|/3sdd - 1| < 7 x 10~ 4 . (2.24) 

This constraint on f^Edd comes only through the combination of «o an d /?o m Eq. ( ^.2l| ) . Therefore the smaller the ao 
is, the looser it constrains /3q. Very recently, Damour and Esposito-Farese [ p6| found that binary-pulsar experiments 
constrain /3q to 

A) > -5, (2.25) 
independently of the value of ao for a particular choice of the functional form a(ip). 



B. Scalar Gravitational Wave 



We show that a scalar mode as well as tensor modes propagates as gravitational waves. We consider linear pertur- 
bations from the Minkowskian metric rj^ v and constant scalar field 0o- 

9nv = iliiv + h liU , (2.26) 
= 00 + 50. (2.27) 



Here afte r in this section we raise and drop tensor indices by the flat metric jf 1 " and ?7 M „. The field equations ( |2.2| ) 
and (2.3) in the linear perturbation are given by |l3| ] 

Sir 1 

= —T^ + —(50.^- v ^n6cj>), (2.28) 

00 00 

□ 50 = 8?r r, (2.29) 
v 3 + 2w V ; 

where cjq = w(0o) an d ft- = h a a . 
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We introduce 



(2.30) 



Then Eqs. (2.28) and (2.29) become 



,<*P_ a P\ 



16tt 



□ r 



— — — T, 

3 + 2wo 



where 9 = 9 a a . We define M „ and h^ v by 



and 



We use the gauge condition defined by 



(2.31) 
(2.32) 

(2.33) 
(2.34) 

(2.35) 



Then the field equations are separated as 



no,. 



16tt 

00 

8?r 



3 + 2w 



-T. 



(2.36) 
(2.37) 



In the absence of matter these equations are wave equations. We can use the remaining degrees of gauge freedom to 
make transverse-traceless. Then the metric perturbations are separated out into the plus mode, cross mode and 
scalar mode. The form of the metric perturbation of the plane wave is expressed by 



'0 0* 

h+ h x 

h x -h+ 

,0 0, 



/-l 0* 

10 

10 

V 1, 



(2.38) 



where h + , h x represent the plus and cross mode, respectively. We can detect the scalar mode by a laser interferometric 
gravitational wave detector as shown by Shibata, Nakao and Nakamura (To). If four detectors are available, we can 
distinguish the plus, cross and the scalar modes and determine the direction of the source in principle [fl()|| . 

Prom Eqs. (|2T0|)-(|2~T^). we can relate the scalar field perturbation 8<f> — 4> ~ 4>o in the Brans-Dicke frame to the 
scalar field perturbation Sip = ip — (po in the Einstein frame in the linear approximation as 



5<j) = — 2ao(l + a )5<p>. 



(2.39) 



III. BASIC EQUATIONS 
A. a?o Expansion 

Hereafter we adopt the Einstein frame. We use the Taylor expansion of a(<p) that characterizes the coupling of the 
scalar field and matter as 



i(tp) = a + (3(<p - <p ) + f3 (2) {ip - ip ) 2 + ■ • ■ . 



(3.1) 
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where ao represents a deviation of the theory from general relativity. We assume that |ao| <C 1. Although it may 
exclude interesting non-perturbative effects, such as a spontaneous scalarization we expand the metric tensor, 

scalar field and energy-momentum tensor of matters in terms of ao, 

T* a p = T^f + a t$ + a§t<§ + O(a 3 ), (3.2) 
9* a p = 9% ] + a h§ + <%h% + O(a 3 ), (3.3) 
<p = tp + a <pW +a 2 ipW +0(al). (3.4) 

From the lowest order of ao i n Eq. ( [2.13] ), we obtain 

G$=»«T$. (3.5) 

( E) ( E) 

This means that and are solutions in general relativity. From the next order of ao, we obtain 

G$=87rt$. (3.6) 

This is also the same as the Einstein equations, so that the metric of the scalar-tensor theory deviates from the general 
relativity by 0(oq). Therefore we can determine the scalar field up to O(ao) by solving the wave equation of the 
scalar field as 

nW<p = -4Tra{tp)T( E l (3.7) 
Back reaction of the scalar field on the space-time can be estimated from O(ag) of Eq. J2.13 ) as 

G% = ftrtS + 2 - \ g ™ g ™>» V $ V $) • (3-8) 

This means that the back reaction can be neglected if relations 

<pM ~ 1, (3.9) 

and 

<4\l%fi\ « (3-io) 

(2) 

are satisfied. Here we have defined the energy-momentum tensor of the scalar field as 

T (2) =lLwJ 1 )_L( B ) fl !% (1 ,W !J (i) , l (-311) 



B. Oppenheimer-Snyder Collapse 



We use the Oppenheimer-Snyder model 27 as the background space-time and matter in which the scalar field 
evolves. The Oppenheimer-Snyder solution describes a gravitational collapse of a homogeneous spherical dust ball. 
In the interior of the dust one can write the line element in the closed Fricdmann form 

ds 2 = -dr 2 + a(T) 2 (d X 2 + sin 2 X dfl 2 ) (3.12) 
= a(n) 2 (-dn 2 + d X 2 + sin 2 X dtt 2 ), (3.13) 
dfl 2 = d9 2 + sin 2 6d<j) 2 , (3.14) 



where 



a (v) = q^C 1 + cos??), (3.15) 
t(v) = 2 a o(^7 + sin 77). (3.16) 
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The density of the dust is given by 



8tt 8nao 



> ! 2 (1 + cosr/) 



The ranges of r\ and x are 



and 



< T) < 7T, 



0<X<Xo < 



(3.17) 



(3.18) 



(3.19) 



Let r s (i) be the circumferential radius of the stellar surface. Then in the exterior of the dust (r > r s (t)), the 
space-time is expressed by the Schwarzschild metric as 



ds 2 = -(l- — ) dt 2 
r 



1 - 



2M 



V 



r J 



dr 2 



r 2 dn 2 . 



(3.20) 



The matching conditions of the interior and exterior solutions are that the proper circumferential radius of stellar 
surface should agree and that the stellar surface should move on a geodesic. These matching conditions are expressed 
as 



r s = 0(77) sin xo, 

1 , 
M = -a sm Xo, 



(3.21) 
(3.22) 



t = 2M In 



£ML _ I] 

v 2M L ) 



tan | 



sin r)) 



(3.23) 



where r s0 = r s (t = 0). 



C. Wave Equation 



We rewrite Eq. (3.7) under the background of the Oppenheimer-Snyder collapse. We use the metric in the form 
(3.14) and (3.20). The wave equation in the interior ( < x < Xo) is given by 



i_d_ f 2 d5tp\ 
a 2 drj \ drj J 



1 d 



dSip\ 
'— — 



Aita(ip)p, 



a'- ^ a" or) \ Or] J sin' x ®X 
while in the exterior (r > r s (t)) it is given by 

\ r J ot £ r z Or y \ r J Or J 
Here in place of Sip we define £ by 

^ _ \ asinxSip (interior) 

1 rSp (exterior). 



Substituting Eq. (3.26) into Eqs. (3.24) and (3.25) , we have 



d 2 ( [ d 2 ( _ 
drj 2 dx 2 
d 2 C d 2 C 2M 



1 



a 
a 



+ 4:Tra(ip) pa 3 sin x (interior), 



m 2 



dr 2 



2M 



£ (exterior), 



(3.24) 

(3.25) 
(3.26) 

(3.27) 
(3.28) 
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where r* is the tortoise coordinate defined by 



(3.29) 



and t he prime stands for differentiation with respect to 77. Using the null coordinates, we can rewrite Eqs. (3.27) and 
dJf ) as 



duS~v = 4 1 + ~ ) C " 8^° SmX (mten0r) ' 



d 2 ( _ M 
dudv 2r 3 

where the null coordinates are defined by 



in the interior and 



2M\ 

1 Q (exterior), 

r 



u = r}-x, 
v = r) + x, 



u = t — r*, 
v = t + r*, 



(3.30) 
(3.31) 



(3.32) 
(3.33) 



(3.34) 
(3.35) 



in the exterior, and we have used Eq. (3.17) to rewrite p. 

The boundary condition at the center of the dust is that the derivative of the scalar field in the radial direction 
should be zero, i.e., 



f^ = at X = 0. 
ox 



(3.36) 



The junction condition on the stellar surface is that ip and its derivative in the direction normal to the boundary 
surface should be continuous, i.e., 



Sf\in 



(3.37) 
(3.38) 



a t X = Xo (interior) and r = r s (t) (exterior), where n M is the normal vector to the boundary surface of the dust. 

For simplicity we take the initial condition that ip = ipo and the time derivative of p> vanishes on the initial 
hypersurface r\ = and t — 0. Namely in the interior of the dust we set 



Sep = 
d8(p 
dr] 



0, 



at 77 = , and in the exterior of the dust 



Sip = 0, 
38 <p 

~~dT 



0. 



(3.39) 
(3.40) 

(3.41) 
(3.42) 



on the first ray u = uq = — r* s (/j = 0), because ip — ipo is the static solution in the exterior region that satisfies the 
boundary condition in the asymptotic region. We can regard this as the initial condition of the situation that for 
t < a highly relativistic star is in equilibrium between the pressure and the gravity, and at t = the pressure is 
switched off and the homogeneous dust begins to collapse, because, for highly relativistic matter, the trace of the 
energy-momentum tensor satisfies 



T = -p + 3P = 0. 



(3.43) 
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IV. NUMERICAL RESULTS 



We divide the Oppenheimer-Snyder space-time to three regions (A), (B) and (C) as shown in Fig. 1 according 
to Cunningham, Price and Moncrief Jl7| . We use null coordinates in the numerical calculation in order to observe 
propagation of the scalar field to the asymptotic region. The details of our numerical calculation will be described in 
Appendix |A|. As for the coupling function, for simplicity, we assume that 

a{v>) = a + P(<p - <Po)- (4-1) 



A. Brans-Dicke Theory 

1. Wave Form in the Wave Zone 

We shall see the behavior of the scalar field in the Brans-Dicke theory (i.e., (3 = 0) to obtain an overall picture of 
the time evolution of the scalar field. Fig. 2 shows the wave form of the scalar gravitational wave by an observer 
at r = 100M in the Brans-Dicke theory from the collapse of the dust whose initial radius r S Q is 10M. The ordinate 
is ( — rdip. The solution is proportional to the parameter cyq and hereafter we shall normalize ( by ao — —0.0316 
corresponding to u> = 500. The first ray reaches the present observer at t ~ 95A/ and then the scalar field increases 
from the initial value ipo to the peak. The amplitude at the peak is estimated as [ flO| 

^~224£. (4.2) 
r 

After the peak, the value of the scalar field ip decreases below the initial constant asymptotic value po and increases 
monotonically toward ipo. A frequency of this last bounce below ipo nearly equals that of the spherically symmetric, 
fundamental quasi- normal mode pi| . This result is in a good agreement with Shibata, Nakao and Nakamura's result 
of full numerical simulations [[l0| . This also shows that the back reaction of the scalar field on the space-time can be 
neglected when lu + 3/2 > 5, i.e., |ao| ^5 0.316. 

2. Time Evolution of the Scalar Field 

Next we shall see the numerical solution £ in the interior of the dust (the region (A)). The result is shown in Figs. 
3(a-c). Fig. 3(a) shows the scalar field in the interior of the dust (the region (A)) in Fig. 1. The initial dust radius 
r s0 is 10M. The case of the initial radius r s0 = 10M is very instructive and therefore we shall first consider this case. 
The abscissas are the null coordinates u = i] — x an( l v = V + X- Fig. 3(b) shows the time evolution of the scalar 
field seen by comoving observers. The abscissa is the conformal time rj and the numbers attached to the curves are 
the values of the fixed radial coordinates x of the comoving observers. Fig. 3(c) shows the time evolution of the 
scalar field configuration on the spacelike hypersurface of the conformal time n = const.. The abscissa is the radial 
coordinate \ an d the numbers attached to the curves are the value of n. 

As already mentioned, at n = the scalar field is (fa, i.e., £ = everywhere. Then ip increases homogeneously in 
the central region (u ~ v) because the source of the scalar field is the homogeneous dust ball and the information of 
the surface of the dust ball has not yet arrived at the central region in the early stage. The information from the 
surface of the dust ball propagates inward at the speed of light and reaches the center at the time n = xo- After the 
reflection at the center, the configuration of the scalar field in the interior of the dust settles to a quasi-static solution 
{dQ/dn ~ 0) along the outgoing null ray u ~ xo and after that the scalar field evolves in a quasi-static manner. 
Finally t he scala r field falls inside the event horizon, keeping its value finite. 



In Sec. IV A3, we will see that the quasi-static configuration is realized only when the initial radius of the dust ball 
is large enough. This fact suggests that the quasi-static evolution appears in the case xo ^ 1: which means that the 
gravity is weak in the initial configuration. In order to confirm this expectation, we shall assume that xo <C 1 and 
that the interior solution £ = £j„ can be expanded with respect to xo- Then we can obtain a consistent quasi-static 
solution up to O(xo), which is given by 

(in = jCi a x(x 2 - 3Xo)> ( 4 - 3 ) 
(ex = -^a a xl ^ -a M. (4.4) 







The derivation of this solution will be described in Appendix |c[ We plot the interior solution (4.3) in Fig. 3(c) with 



the label "QS". We see that the approximate solution agrees with the numerical solutions very well for the case of 
(3 = even when the dust surface is rather close to the event horizon. When the surface of the dust ball is close to 
the event horizon, the amplitude of the scalar field in the interior region is 

6^-^, (4.5) 
r s 

and the ratio of the energy density of the scalar field to that of the dust remains 0(a\). 

The numerical solution C, = rSip in the exterior of the dust is shown in Figs. 4 (the region (B)) and 5 (the region 
(C)). The set of null coordinates u , v of Fig. 4 is a continuous extension of the interior characteristic coordinates 
obtained by relabeling the rays u = const, and v = const.. The details of this coordinate extension will be described 
in Appendix |b| In the region (B) the scalar field increases and settles to the quasi-static configuration on the null ray 
u ~ xo- hi the region (C) the outgoing radiation of the scalar field can be seen. If we retrieve G and c, the units of 
time and length become 

^=4.93xl0-#- sec, (4.6) 



c 3 M 

Vrr 



GM M , 

— = 1.48— km, (4.7) 



respectively. 

As seen in Figs. 4 and 5, first the scalar field increases from the initial value ipo due to the presence of the dust. 
When the event horizon is formed, the field inside cannot affect the field outside. At that moment no wave propagates 
outwardly from the region (B) to the region (C). In the wave zone the scalar field approaches the asymptotic value 
tpo after the wave has passed the observer at r = const.. This is consistent with the Hawking's theorem which states 
that the black hole has no scalar charge jl5| . In Fig. 5 a steep cliff can be seen on the boundary v = vq between the 
regions (B) and (C) but this is only due to the coordinate singularity of the Schwarzschild coordinates and does not 
mean any singular behavior of the scalar field. 



3. Dependence on the Initial Radius of the Dust 

In Fig. 6, wave forms of the scalar gravitational waves from the collapse of the dust whose initial radii r S Q are 
3M, AM, 6M, 8M, 10M and 20M in the Brans-Dicke theory are plotted. As we can see in Fig. 6, peaks are nearly 
symmetric for r s0 < AM but not for r s0 > AM. Fig. 7 shows the relation between the peak amplitude of the wave 
form and the initial dust radius. As seen from Fig. 7 , the amplitudes at the peaks are nearly |cto | Af/ r for r s o ^ AM 
but smaller than \ao\M/r for r S Q < AM. These can be explained as follows. The information on the stellar surface 
reaches the center of the dust first at r\ = xo- On the other hand the information at the center at n > ir — 3xo cannot 
travel to the exterior because it is inside the event horizon. Therefore for the dust of xo ^ 7r /4, i.e., r s o < AM, the 
scalar field cannot reach the quasi-static evolution until the event horizon is formed. So a mild downhill after the 
peak, which is an evidence of the quasi-static evolution in the interior of the dust , cannot be seen for r S Q < AM. On 
the other hand, for r s o ^ AM, a mild down hill after the peak implies the quasi-static evolution in the interior of the 
dust. 

For r s o > AM, since the initial ingoing null ray from the stellar surface bounces back to the surface rj ~ 2xo, the 
time At pea k that it takes for the scalar field to reach the peak from its rise is estimated as 



At peak ~ t(n = 2xo) 

cot xo + tan xo 



2M In 



cot xo - tan xo 



2M cot xo 



2xo 



1 



2 sin xo 



(2xo + sin 2xo ) 



(4.8) 



This shows an agreement with the numerical results in Fig. 6 within an accuracy of ~ 10 %. For example, using Eq. 
(4.8), we can estimate Ai pea fe to be ~ 40M for r s0 = 20M. 

Fig. 8 shows the relation between an initial dust radius r s o and a characteristic frequency f c which is defined as the 
frequency at which the energy spectrum of the wave form takes the maximal value. For r s o ^ 8M, the characteristic 
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frequency f c is approximately in proportion to the inverse of the free-fall time of the initial dust, which is shown by 
Shibata, Nakao and Nakamura |10| as 

f c cx J- ~ JpS. (4.9) 

For r s0 < 8M, though there is no such simple relation because of the effect of the space-time curvature, f c would 
be determined by the initial radius and mass of the dust. The last bounce below the asymptotic value <po reflects a 
quasi-normal oscillation. Fig. 9 shows the relation between the initial dust radius and the time scale of the bounce, 
which is defined as the half value width of that bounce. As shown in Fig. 9, the time scale of this bounce is almost 
independent from the initial dust radius and nearly equals to the period of the quasi-normal oscillation f28[| . 



B. Dependence of the results on the Parameter in the Theory 

In Figs. 10-12 the wave forms of the scalar gravitational waves seen by an observer at r = 100M are shown, for 
various (3, from the collapse of the dust whose initial radius is fixed to 10M. In each figure, the scalar field approaches 
the asymptotic value after the peak, the quasi-normal mode of the scalar wave having passed the point of the observer, 
which is consistent with the black hole no hair conjecture p9[ . 

First we consider the case of positive (3. The amplitude of the scalar gravitational wave decreases as (3 increases. 
As seen from Figs. 10 and 11 there appears another oscillation mode whose frequency is characteristic of the relevant 
gravitational theory and the dust density. The larger /3 is, the shorter is the period of this oscillation mode and the 
more the scalar field oscillates before it approaches the asymptotic value. This can be also seen from Fig. 13, which 
shows the time of the first minimum from the reach of the first ray. This oscillation mode reflects an oscillation mode 
in the interior of the dust seen in Fig. 14, which shows the time evolution of the scalar field observed by comoving 
observers in the region(A) for (3 = 50 . The ratio of the energy density of the scalar field to that of the dust remains 
within ~ O(oq). So, if |ao| -C 1, the back reaction of the scalar field on the space-time is negligible. 

Second we consider the case of negative j3. The amplitude of the scalar wave increases as \(3\ increases, as seen in 
Fig. 12. This reflects an exponential growth of the scalar field in the interior of the dust seen in Fig. 15, which shows 
the time evolution of the scalar field observed by comoving observers in the region(A) for (3 = —5 . The ratio of the 
energy of the scalar field to that of the dust is ~ 10 x 0(oq) for /? = — 1, ~ 100 x 0{oq) for (3 = — 5 and ~ 10 4 x 0(oq) 
for (3 = —10. Therefore, if we consider a finite value of ao, there is a critical value of (3 for which the back reaction of 
the scalar field cannot be neglected. 



As we have seen in Sec. [IA, in fact, the solar-system test experiments for the post-Newtonian order deviation 



from general relativity do not constrain directly the value of (3 because (3 is constrained only through the combination 



with ao in Eq. (2.21) (however, see |26[). So the wave form for (3, the value of which is allowed by the solar-system 



experiments, may be quite different from that in the Brans-Dicke theory. 



V. DISCUSSIONS AND SUMMARY 

In the Brans-Dicke theory the peak amplitude of the scalar gravitational wave is given by 

500 \ / M \ /10Mpc\ 



ty-oo^-lO-" ( — ) ( — 1 (— rH, (5-1) 



and the characteristic frequency is given by 



/c^xlosf^y'V^r^V (5.2) 



M Q J \15km/ 

if r s0 ^ 4M. For this frequency the sensitivity of the first LIGO would be h ~ 10~ 21 and that of the advanced LIGO 
would be h ~ 10 -22 . If oj ~ 500, the advanced LIGO may detect the scalar gravitational wave from collapse of a nearly 



spherical mass ~ 1M & and initial radius ~ 15km at the distance ~ IMpc from us. If the delayed collapse []3Q 31 
after a supernova would happen in our Galaxy, we can test the Brans-Dicke theory by the advanced LIGO up to 
lo < 10 5 . Even for a neutron star formation, scalar gravitational wave may be emitted because the progenitor would 
lose a considerable part of its scalar mass in the process. The amplitude of the wave would tell us the information 
of its self-gravitational energy. If a wave form of the scalar gravitational wave from gravitational collapse is obtained 
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observationally, we can determine its amplitude, characteristic frequency and quasi-normal mode frequency. We can 
then determine the source mass because the frequency of the quasi-normal mode is in proportion to the inverse of 
the mass. If we know, in addition, the distance to the source by another method, we can determine the Brans-Dicke 
parameter lu. Moreover, we can determine the initial dust radius from its characteristic frequency. In the case of 
delayed collapse after the supernova, the core radius would correspond to it. The information of the initial radius 
would constrain the possible high-density equation of state. 

Next we consider the parameter dependence of the wave form in the scalar-tensor theory. The wave form greatly 
depends on 0. If the space-time is flat and the density is homogeneous, a homogeneous solution which describes 
harmonic oscillations at the period \Jit/(l3p) for j3 > and an exponential increase at the e-fold time l/(2y / n\f]\p) 
for (3 < appears. Such a mode would appear in the wave form of the scalar gravitational wave. For (3 > the back 
reaction of the scalar field on the space-time is negligible, and so we can safely expand the theory around general 
relativity. Therefore we can determine the wave form of the scalar gravitational wave within our approximation. For 
(3 < the back reaction is not negligible and the full calculation that deals with the dynamics of the metric, matter 



and scalar field is needed in order to predict the correct wave form ||32fl . We also see in Eq. (3/7) that the scalar field 
gets an effective mass m 2 — ~Anj3T. For (3 > 0, the scalar field can be regarded as being effectively massive since 
T = —p + 3P < , and for (3 < 0, the scalar field gets effectively a tachyonic mass. 

Finally we make a crude estimate of the condition in which the scalar mode of gravitational waves radiated from 
its collapse dominates over the tensor mode. We consider collapse of an oblate spheroid of mass M. Then the tensor 
mode of the wave at distance r takes the form 

h T ~-Q ~ -M(a 2 - c 2 )-, (5.3) 
r r 

where a and c are the semi -major axis and semi-minor axis respectively. While the scalar mode is 

*.«±^, (5.4) 
U) r 

where Q is the reduced quadrupole moment and the dot means a time derivative. If we replace the time derivative 
with the free-fall time tff , then the ratio of hr to hs is 

, 2cjM(a 2 ~c 2 ) , 2 2 M ,/w\/10M\ e 2 ., c , 

hT,hs ~ ~m t^~f ~ 1Ma " c } ^ ~ 10 (500) (— ) 7T=W (5 - 5) 

For example, for u> — 500 and a = 10M the scalar mode dominates over the tensor mode if the eccentricity e < 0.03. 
Thus the scalar mode would dominate in almost spherical collapse. 

We have examined a scalar gravitational wave from a spherically symmetric and homogeneous dust collapse in the 
Brans-Dicke theory and scalar-tensor theories of gravity, where we have neglected the back reaction of the scalar field 
on the space-time. 

(1) In the Brans-Dicke theory the back reaction of the scalar field on the space-time remains 0(l/ui). So if w ^ 1, 
it is negligible. 

(2) In the Brans-Dicke theory the amplitude of the scalar gravitational wave is given by ~ M/(u>r) and the 
characteristic frequency is given by ~ 3 x lO 3 (Af/M0) 1 / 2 (r s o/15km) _3 / 2 Hz. The characteristic frequency of the 
scalar gravitational wave is in proportion to the inverse of the free-fall time of the initial dust sphere. The last bounce 
reflects the fundamental quasi-normal mode of the black hole. If r s g > AM, the scalar field in the interior of the dust 
reaches a quasi-static configuration and thereafter evolves quasi-statically. Therefore the wave form after the peak 
goes down slowly. If r s o < AM, the scalar field in the interior of the dust does not reach the quasi-static configuration 
and so the wave form after the peak goes down more quickly. The peak amplitude of the scalar gravitational wave is 
smaller than Mj (tor) in this case. 



(3) Both in the Brans-Dicke theory and in the scalar-tensor theory given by Eq. (4.1), the scalar field approaches its 
asymptotic value after the scalar wave radiation has passed the observer at r = const.. In our numerical calculation 
we do not find any singularity in the scalar field at the event horizon. 

(4) In the Brans-Dicke theory we can determine a mass, initial radius of the dust and the Brans-Dicke parameter 
from the wave form of the scalar gravitational wave and the distance to its source. 

(5) We calculate the scalar gravitational wave from the dust collapse in the scalar-tensor theory given by Eq. (p~l|). 
The wave form greatly depends on the parameter (3. If j3 > 0, the amplitude of the scalar gravitational wave is 
suppressed and a new oscillation mode whose time scale is ~ l/\/j3~p appears. The back reaction of the scalar field 
on the space-time is within ~ 0{ocq). If (3 < 0, the amplitude of the scalar gravitational wave is enhanced and its 
back reaction would not be negligible. Hence we have to solve the whole field equations fully numerically, which is 
the subject of our next paper |32|. 
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APPENDIX A: NUMERICAL CODE 



The wave equations we calculate take the form. 

dq 
du 

21 
du 
dp 

dv 

at 



:A(u,v,£,p,q), (Al) 

P, (A2) 

: A( u > v >^p>q)> ( A3 ) 

: q, (A4) 



Our numerical code is similar to that of Hamade and Stewart |33|], but we slightly modify their finite differential 
equations. We solve simultaneous partial differential equations given by 

^ = G(y,z). (A6) 

We determine y n , z n at the point n(u, v) when y w and z w are given at the point w(w, v — h) and y e , z e at the point 
e(it — h,v). At the first step we calculate 

lln = Ve + hF(y e , z e ) (A7) 

z n = z w + -h(G(y w ,z w ) + G(y ni z n )), (A8) 

and at the second step 

Vn = tj(2M + Ve + hF(y n , z n )), (A9) 
z n = ^ j^n + z w + ^h(G(y w ,z w ) + G(y n , z n ))X . (A10) 
In general this is an implicit scheme, but since the rig ht ha n d side of the wave equation we have calculated does not 



contain any first derivative of the field as seen Eqs. ( |3.3C| ), ( 3.31 ), it becomes an explicit scheme. The error of the 
one set of these two steps is estimated as 0(h 3 ). 

This code was checked through the following tests. Wc calculated a spherically symmetric, ingoing wave in flat 
space-time which bounces at the center, and compared it with an analytic solution. The error of this test is within 
~ 0.5%. We calculated scalar gravitational waves from the Oppenheimer-Snyder collapse in the Brans-Dicke theory 
and compared them to the results for u> = 500 of the full simulation by Shibata, Nakao and Nakamura |l(| . The two 
results are agreed within an accuracy of several percent. These tests make us sure that our numerical code works 
well. It took about two hours to compute one model on a Sun SPARC Station 20 (50MHz). 



APPENDIX B: COORDINATE EXTENSION 



It is convenient for our numerical calculation to divide the exterior region into the two regions (B) and (C) in Fig. 1 
by the inward light ray v = vq which reaches the stellar surface on the event horizon. We extend the interior comoving 
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coordinates to the exterior region following Cunningham, Price and Moncrief |L7]]. We relabel light rays u = const. 
and v = const, by the interior coordinates u, v at which these light rays reach the stellar surface. This coordinate 
transformation is expressed explicitly as 



u = t s (u + xo) - r* s {u + Xo), 
v = t s (v - xo) + r* s (v - xo), 
cot xo + tan 3 



t s {n) = 2M\n 



COt XO ~ t an 2 



+2M cot xo 



1 



r* a {n) = a(n) sinxo + 2Mln 



2 sin 2 xo 
a{n) sinxo 



2M 



(?7 + sin n) 
1 



(Bl) 
(B2) 



(B3) 
(B4) 



This new coordinates are well-behaved everywhere in the space-time except for the central singularity of the space- 
time. 

In these coordinates, the junction conditions on the surface become particularly simple because the vector normal 
to the surface is given 



n 



1 f d_ _ d_ 

a \dv du 



(B5) 



APPENDIX C: QUASI-STATIC SOLUTION 

We assume that Xo ^ 1 an d that the interior solution £ = Q n can be expanded with respect to xo as 

oo 

Cin = E4 n) . (CI) 



71 = 1 



where Q"' = O(xo)- Since x is small, the interior solution Q n will be well approximated by the truncated Taylor 
series of x an d therefore we should impose 



(n) 



d x r > 



o{ X r m )- 



(C2) 



Under the assumption of the quasi-static evolution phase, d(i n /dn = d 2 (i n /dn 2 = 0, the wave equation ( 3.27 ), in 
the case of /3 = 0, is written up to O(xo) as 



(C3) 
(C4) 
(C5) 
(C6) 
(C7) 
(C8) 



d 2 ( (0) 

dx 2 


= 0, 




5 2 C (1) 
dx 2 


= 0, 




d 2 c {2) 

dx 2 


ao JO) 
2a W ' 




d 2 c {3) 

dx 2 


a-oJi) 
2a W 


3 

- 2 a oa X 5 


d 2 c {4) 

dx 2 


_ ao J2) 
2a W ' 




d 2 c {5) 

dx 2 


ao J3) 
2a 


1 3 
- ^aoaox 



where we have used Eqs. ( 3.15 ) and ( 3.17 ) to rewrite a" and p 
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On the other hand, we consider only the nontrivial, lowest-order solution for the exterior region. Since 2M/r s - 
(ao/a) sin 2 xo = O(xo)i as long as ao/a is not so large, the lowest order wave equation for the exterior solution £ = ( e 
becomes 



d Cex . d Cea 

dt 2 dr 2 



0. 



(C9) 



Under the assumption of quasi-static evolution, the matching condition between the interior and the exterior solutions 
at the dust surface r — r s , that is, x — Xo is given by 



Cea 
d(ex 

dt 

d(ex 

dr 



Cin > 
1 . 



a 
1 

a 



■smxo 



a 



I , a ■ 2 

- 1 sin xo 

* a 



9x 
-l 



1 

a 

cos xo- 



no 
a 



1 



d(in 

dx 



dx a dx 



Xo - 
[1 



0(xl)}, 
-O(Xo)]- 



From the above equation, we find that 



\d( ex /dt\ 
\d( ex /dr\ 



O(xo) 



at 



(CIO) 
(Cll) 

(C12) 
(C13) 



This means that th e ex terior solution satisfies dC, ex /dt — in the lowest order of xo as long as ao/a <C 1/xo and 
therefore from Eq. (C9), we obtain 



where we have imposed the boundary condition 
condition ( p!2 ) in the lowest order becomes 

dQ, 



Ccx ^ex COTiSt. , 

+ const, for r 



(C14) 

co. From the above equation, the matching 



dx 







at 



X = Xo- 



Now we return to the interior solution. The solutions for (n = 0, 1) are given by 

and 6n } =C^ + D^ X , 



jo) = r (o) 

>in in 



(C15) 



(C16) 



where (n — 0, 1) and are the integration constants and the order of those are, respectively, d|„ = 0(1), 



C 



(i) _ 

in 



O(xo)- However from the regularity at the origin, x = 0, i.e., Q n oc x f°r X ~~ * 0, and vanish, while 



.(i) 



is determined by the matching condition ( |C15| ) and we find also vanishes. As a result, we get 



= C 



Using the above results, we obtain the solutions of the equation for (n = 2, 3) as 



c, 



(2) 



0. 



and 



c 



(3) 



1 



r^oa-oX" 



D (3) v 



(C17) 



(C18) 



(3) 

where we have used the reg ulari ty condition at the origin and D ir / is the integration constant which is determined 
by the matching condition ( C15 ). Hence th e no n-trivial solution for the interior region appear s at O(xo) an d using 
this solution, C ex is determined from Eq. (C1C). Then we obtain the solution given by Eqs. ( |4.3| ) and (4.4). This 
solution is clearly time independent and is consistent with the assumption of quasi-static evolution. From Eq. ( |C7| ) , 

— 0, while Eq. (C8) for (jj^ becomes 

(C19) 



the first equation of Eq. (|C18|) and the regularity at the origin, we obtain Q, 



2^(5) 
in 



dx 2 



a 2 1 
-^<*ox(x 2 ~ 3xo) - -a aoX 3 - 



The above equation shows that Q n depends on time r\ through the first term on the right-hand side and hence the 
assumption of the quasi-s tatic e volu tion is no longer consistent with the solution above O(xo)- ^ n t ne case of j3 ^ 0, 
the quasi-static solution (4.3), (4.4) is also valid up to O(xo)- Here we should again note that our approximation 
based on the smallness of xo is valid only when ao/a is not so large, i.e., ao/a <C 1/xo and further ao/{2a) <C 1/xo- 
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FIGURE CAPTION 



Fig.l. Regions of the Oppcnhcimcr-Snyder space-time for constant 9 and </>, expressed in characteristic coordinates. 
In the "stationary region" the initial stationary field configuration remains a stationary solution. 

Fig. 2. The wave form of the scalar gravitational wave in the Brans-Dicke theory at r = 100M. The initial dust radius 
r s o is 10M. The ordinate is ( — rSip. The abscissa is the time t from the beginning of the collapse at t = 0. 

Fig. 3. The scalar field in the interior of the dust (the region (A) of Fig. 1) in the Brans-Dicke theory The initial 
dust radius r s0 is 10M. (a) The ordinate is £ = asinxSip. The abscissas are the null coordinates u = n — \ 
and v = rj + x- (b) The time evolution of the scalar field seen by comoving observers. The numbers attached to 
the curves are the radial coordinates \ °f the comoving observers. The ordinate is £ = a sin x an d the abscissa 
is the conformal time r\. (c) The scalar field configuration on the spacelike hypersurface of the conformal time 
rj = const.. The abscissa is the radial coordinate x- The numbers attached to the curves are the value of r\. 
The curve with the label "QS" shows the approximate quasi-static solution obtained in the text. See the text 
for further details. 

Fig. 4. The scalar field in the exterior of the dust (the region (B) of Fig. 1) in the Brans-Dicke theory. The initial 
dust radius r s0 is 10M. The ordinate is ( = rSip. The abscissas are the null coordinates u and v extended from 
the interior region to the exterior. 

Fig. 5. The scalar field in the exterior of the dust (the region (C) of Fig. 1) in the Brans-Dicke theory. The initial dust 
radius r s0 is 10M. The ordinate is C — rSip. The abscissas are the null coordinates u = t — r* and v = t + r*. 

Fig. 6. The wave forms of the scalar gravitational waves from the collapses of various initial radii in the Brans-Dicke 
theory. The ordinate and abscissa are the same as Fig. 2. The numbers attached the curves represent the value 
of the initial radii respectively. 

Fig. 7. Relation between the peak amplitude and initial dust radius r s o in the Brans-Dicke theory. 

Fig. 8. Relation between the characteristic frequency f c and initial dust radius r s o in the Brans-Dicke theory. 

Fig. 9. Relation between the time scale of the last bounce and initial dust radius r s0 . The ordinate is the time scale 
of the last oscillation which is defined as the half-value width of that bounce. 

Fig. 10. The wave form of the scalar gravitational wave for various (3 at r = 100M. The numbers attached to the 
curves represent the value of (3 respectively. The initial dust radius r s0 is fixed to 10M. 

Fig. 11. Same as Fig. 10, but for different [3. 

Fig. 12. Same as Fig. 10, but for different (3. 

Fig. 13. The time of the first minimum of the wave form from the reach of the first ray as a function of (3 > 0. 

Fig. 14. The time evolution of the scalar field seen by comoving observers in the interior of the dust (the region (A) 
of Fig. 1) for (3 = 50 . The ordinate and abscissa are same as Fig. 3(b). 

Fig. 15. The time evolution of the scalar field seen by comving observers in the interior of the dust (the region (A) of 
Fig. 1) for (3 = — 5 . The ordinate and abscissa are same as Fig. 3(b). 
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